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ABSTRACT 

The role of a photospheric component and of pair breakdown is examined in 
the internal shock model of gamma-ray bursts. We discuss some of the mechanisms 
by which they would produce anomalously steep low energy slopes, X-ray excesses 
and preferred energy breaks. Sub-relativistic comptonization should dominate in 
high comoving luminosity bursts with high baryon load, while synchrotron radiation 
dominates the power law component in bursts which have lower comoving luminosity 
or have moderate to low baryon loads. A photosphere leading to steep low energy 
spectral slopes should be prominent in the lowest baryon load cases. 
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1. Introduction 

The standard fireball shock scenario for gamma ray bursts (GRB) assumes a synchrotron 
and/or an inverse Compton (IC) spectrum, in good general agreement with a number of 
observations (Tavani 1996; Cohen et al. 1997, Panaitescu, Spada & Meszaros 1999). The physical 
motivation for this scenario is strong; it would indeed be surprising if the expansion of the ejecta 
from the huge energies inferred in GRB did not lead to shocks, where synchrotron and IC radiation 
play a significant role. One need only remind oneself of the dominant role of these effects in AGN 
jets and supernova remnants. However, in a significant fraction of bursts there is evidence in 
the 1-10 keV range for spectral intensity slopes steeper than 1/3 (photon number slopes flatter 
than -2/3), as well as in some cases for a soft X-ray excess over the extrapolated power law 
from higher energies (Preece et aZ.1998, 1999; Crider et a/.1997, etc.), and this has served as 
the motivation for considering a thermal (e.g. Liang (1997) or a thermal/nonthermal (Liang et 
al. 1999) comptonization mechanism. While an astrophysical model where this mechanism would 
arise naturally has been left largely unspecified, Ghisellini & Celotti 1999 have pointed out that 
if internal shocks lead to high compactness parameters and pair formation, this could naturally 
produce conditions where the pair temperature and the scattering opacity are self-regulated at 
values favoring comptonization. 

There is also a trend in current analyses of GRB spectra indicating that the apparent 
clustering of spectral break energies in the 50 keV-1 MeV range is not due to observational 
selection effects (e.g. Preece et aZ.1998; Brainerd et a/.1998; see however Dermer et a/.1999a). 
Models to explain a preferred break physically, e.g. through a Compton attenuation model 
(Brainerd et al. 1998) require reprocessing by an external medium whose column density adjusts 
itself near a few g cm -2 . More recently a preferred break has been attributed to the blackbody 
peak of the fireball photosphere when this occurs at the comoving pair recombination temperature 
in the accelerating regime, which is blueshifted to the appropriate observer frame energy (Eichler 
& Levinson 1999), and in this case the Rayleigh- Jeans portion of the photosphere provides a steep 
low energy spectral slope. To get a photosphere as close in as to occur at the pair annihilation 
temperature requires an extremely low baryon load outflow; the presence of a high energy power 
law extending to GeV then requires a separate explanation. 

At the other extreme of large baryon load outflows, Thompson (1994) has considered scattering 
off MHD turbulence in photospheres in the coasting regime, which boosts the adiabatically cooled 
thermal photons near the photospheric peak up to a larger break energy, leading to canonical 
energy slopes of and -1 (photon slopes -1 and -2) below and above the break. An alternative 
scenario is that of Ghisellini & Celotti (1999), who invoke pair breakdown in high compactness 
shocks above the photosphere to achieve a self-regulating nonrelativistic lepton temperature 
and scattering optical depth of order a few, leading to a thermal comptonization spectrum 
with standardized features. While this does not explain steep low energy slopes it does provide 
nonthermal spectral energy power laws of slope ~ close to those of synchrotron, as well as a 
standard temperature which would lead to a preferred break in the observer frame, if the bulk 



-3- 



Lorentz factor is in a narrow range. 

In this paper we synthesize and extend some of these ideas within the framework of the 
standard fireball internal shock model. While in our previous work (e.g. Meszaros , Laguna & 
Rees 1993, Rees & Meszaros 1994) we considered photospheres and pair formation, their thermal 
character, the strong uncompensated photosphere redshift in the coasting phase, and the lack 
of a straightforward way to get from these a power law extending to GeV energies served as 
compelling arguments for concentrating on the synchrotron and inverse Compton mechanisms. 
Here we re-examine the relative roles of photospheres and shocks, as well as those of synchrotron, 
pair breakdown, scattering on MHD waves and comptonization. This leads to a unified picture 
where both shocks and/or a photosphere with a nonthermal component can provide most of the 
luminosity, and where the synchrotron and IC mechanisms in shocks provide the primary spectrum 
or (in high comoving luminosity cases) the trigger for pair breakdown leading to comptonization. 
We investigate the range of burst model parameters over which different mechanisms come to 
the fore, and discuss their role in providing flatter or steeper spectral slopes as well as preferred 
spectral break energies. 



2. Fireball Photospheric Luminosity 

We assume a fireball wind of total luminosity output L Q = 10 52 Ls2erg s _1 expanding from 
some initial radius r Q , which for the sake of argument is normalized to the last stable orbit r Q 
at three Schwarzschild radii around a non-rotating black hole of mass M^h = IO/ziM©, with a 
corresponding Kepler rotation timescale t Q , 

r Q = 6GM bh /c 2 = 0.9 x 10 Vi cm , 

t = 2itrl /2 (2GM bh y 1 / 2 = 3.25 x lO^i s . (1) 

The initial blackbody temperature in electron rest mass units at that radius is 

e o = {k/m e c 2 ){L /^r 2 cT 2 a) l l A = 2.42L^^ 1/2 F- V 2 , (2) 

(or ~ 1.2 MeV) for a 10/ii solar mass BH and an initial bulk Lorentz factor r o > 1 at r Q . As the 
optically thick (adiabatic) wind expands with comoving internal energy e' oc n' 4 / 3 , where n' is 
the comoving baryon number density, the baryon bulk Lorentz factor increases as T oc r and the 
comoving temperature drops ©' oc r _1 . The pairs drop out of equilibrium (Paczyhski 1986, 
1990) at @' p ~ 0.03 (~ 17 keV), at a radius r p where the bulk Lorentz factor has grown linearly to 
a value F p , 

? = ^ = ^ = 7 x lO^V'V' 2 • (3) 

This is the radius of an pair photosphere, above which the scattering optical depth is less than 
unity, unless the wind carries enough baryons to provide an electron scattering photosphere above 
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r p . For a wind baryon load M parameterized by a dimensionless entropy r\ = L/Mc 2 , the baryonic 
electrons lead to a photosphere larger than equation (||) if 77 < rj p (equation |7j ) . As long as the 
wind remains optically thick it is radiation dominated and continues to expand as a relativistic 
gas with r oc r (e.g. Shemi & Piran 1990). Clearly T cannot exceed rj = L a /Mc 2 , and for large 
baryon loads or moderate rj this occurs at a saturation radius r s /r Q = rj/T (for low loads or large 
r\ see however equation (Tl]]). Above the saturation radius, the flow continues to coast with T = 
constant equal to the final value achieved at r s . 

An electron scattering photosphere is defined by t' s = n'YaTr p h/T = 1, where 
n! = (L/A7rr 2 m p c 3 rr]) is the comoving baryon density, Y is the number of electrons per 
baryon and r/T is a typical comoving length. For relatively low values of rj < i]* (defined in 
equation ||) the flow remains optically thick above the saturation radius r s and the photosphere 
arises in the coasting regime T = rj = constant, at a radius r p ^ > r s , (Rees & Meszaros , 1994, 
Thompson 1994), 

' o 

= r-V^A?*) -3 ■ (4) 



Here ??* is the critical value at which 

La T YT 



47rm„c 3 r 



~ 10 3 (L^YTo) 1 ^ , (5) 



which is the wind equivalent of the critical rj of an impulsive fireball photosphere discussed in 
Meszaros , Laguna & Rees 1993, labeled there rj m . 

For low baryon loads where r\ > 7/* a baryonic electron photosphere appears in the accelerating 
portion T oc r of the flow at r p h < r s , 

r ph _ ( La T Y \ X / 3 _ o, v in3fr„Vr-2U/3„-l/3 -1/3 



f o = (4^%^) ' =2.35xl 3(L 52 yr-)V3^ 



r-V(r7/V)~ 1/3 • (6) 

These photospheric radii are shown in Figure 1. For sufficiently high rj, the baryon photospheric 
radius given by equation (|6|) can formally become smaller than the pair photosphere radius of 
equation (||), in which case the latter should be used. The minimum possible photospheric radius 
is therefore achieved at r = r p given by equation (^), requiring extremely large values of rj > rj. p 
(> V*), 

Vp = ix l^L^^YT- 1 ' 2 = 3.75 x lO 3 ^!^ 1 ) 3 / 4 , (7) 

which implies extremely low baryon loads, M < 1.5 x 1O~ 8 M0S -1 . 

The Lorentz factor attained at the photosphere, T p h, grows linearly with rj for 1 < rj < r/* 
(assuming r o = 1), then it decays as r/ -1 / 3 up to rj p and remains constant above that. The 
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comoving temperature Q' ph oc n' 1 / 3 oc r 1 for a photosphere at r p h < r s , while &' p oc r 2 / 3 for 
r ph > r s- The observer-frame photospheric temperature Q p h = Q ph T p h is then 

Qpfe = | (r p h/r s )~ 2/3 = (r?/r?*) 8/3 , for 77 < 77*, r p/l > r s ; 
©o \ 1, for ?? > 77*, rph < r s . 

The observed photospheric thermal luminosity is L p h t oc r 2 r 2 0^ oc r° for r < r s and Lp/^ oc r~ 2 / 3 
for r > r s , hence 

-^pM _ I (r ph /r s y 2 / 3 = (t?/??*) 8 / 3 , for 77 < 7/*, r p h > ?~ s ; ^ 
L \ 1, for 77 > 7/*, r p?i < r s . 



3. Kinetic and Internal Shock Luminosity 

For typical values of 77 £ 10 3 it is clear that the terminal baryon Lorentz factor is F = 77, since 
the photosphere occurs beyond the saturation radius after the baryons are already coasting with 77. 
However the terminal baryon Lorentz factor is less obvious in cases where 77 > 77*. For such values, 
a photosphere occurs in the regime where T oc r, so r p h < r s , but the question is what happens 
to the baryons above this photosphere, and what is the appropriate value of r s . One possibility is 
that the outflow has magnetic fields strong enough that Poynting stresses continue to accelerate 
baryons outside the photosphere. If radiation provides the dominant relativistic pressure, to 
achieve a saturation radius at the value r s /r Q = rj/T would require the baryons to be coupled to 
the radiation out to that radius, beyond the photosphere. Alternatively, it is sometimes assumed 
that the baryons decouple at the photosphere r p h, and coast thereafter with T = T p h = r o (r p h/r ). 
However, the fact that the outflow has become optically thin to scattering means that most 
photons no longer scatter. Nonetheless, most of the electrons above the photosphere can still 
scatter with a decreasing fraction of free-streaming photons, as long as the comoving Compton 
drag time t' dr = m p c 2 /caTu' r is less than the comoving expansion time t' ex = r/cT. The ratio of 
these two times, 

(4/4) = (47rm p c 3 rT 3 /La T ) = (7?,/r o )- 4 (r/r ) 4 , (10) 

exceeds unity above a radius r*/r a = 7/*/r o . Thus for rj > 77* the appropriate saturation Lorentz 
factor is 77* (instead of the larger 77) and the saturation radius is r s /r Q = r*/r Q = 77* /Y a < t]/T . 
Thus, in general the terminal bulk Lorentz factor and the saturation radius are 

r s = min [77, 77*] , (r s /r ) = min [77, 77*]/r , (11) 

where the critical value 77* is given by equation (||). 

The observer-frame kinetic (matter) luminosity of the outflow oc r 2 r 2 n' (kT' + m p c 2 ) oc r 
for r < r s and L/% oc r° for r > r s . For 77 < 77* it is clear that reaches the level L a = Mc 2 rj 
since the terminal bulk Lorentz factor saturates at the initial dimensionless entropy 77. However, 
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for 77 > 77* the terminal Lk can only reach a lower level, since the bulk Lorentz factor saturates at 
the lower value 77* < 77. The terminal value of L/% above r s is then 

Lk = (l, for 77 < 77*, r > r s ; 

£ I (W?7) < 1 > for 77 > 77*, r > tv 

Internal shocks can occur when the flow has variations in the initial 77 or L with consequent 
variations in the terminal Lorentz factors, so shells with different T s catch up with each other. The 
shocks cannot occur at r < r s since in this region both shells accelerate at the same rate T oc r and 
do not catch up. After V has saturated at r > r s , shocks develop at radii ~ 2ct v TiT2 ~ 2ct v T 2 for 
shells whose terminal Lorentz factors differ by Ar = T2 — T\ ~ T as a result of initial variations in 
77 or L over timescales t v = £ v t a < t w , where £ v > 1 and t Q is the minimum dynamic timescale in 
equation ([[]). 

(r sh /r ) = 2.17 x 10%mTl = 2.17 x lO^fT/^) 2 . (13) 

(The factor 21.7 = 2 x 2vr x V3 comes from the factor 2ct v in the shock radius definition, the 
2tt from taking the rotation time at r D rather than crossing time, and y/3 because r Q is at three 
Schwarzschild radii.) 

To produce non-thermal radiation, shocks must occur in an optically thin region, which 
requires 

V > Vsh,m = 1-42 x 10 2 L^ 5 /x- 1/5 ^ 1/5 C V5 . for r ah > r> h , (14) 

in order for a shock to occur above a photosphere which is in the coasting region, r s h > r^ h > r s . 
For £ v = 10 3 corresponding to variability timescales 10 3 i o ~ Is, this can be as low ELS Tjgfi m ~ 35. 
For even higher loads such that 77 < rj s h^ m , the photosphere is further out on the coasting 
regime, and any shocks would occur inside the photosphere. Internal shock radii are shown in 
Figure 1 as a function of 77 for various multiples £ v = t v /t of the Kepler time at the last stable 
orbit t Q . However, from causality considerations shocks may occur at even lower radii, formally 
corresponding to > 1/21.7, since as soon as the r > r s coasting regime is reached shells of 
different 77 can catch up with each other. This would allow smaller regions where the variability 
timescale is as small as r Q /c. 

For very low baryon loads or very high 77 > 77*, the photosphere arises in the accelerating region, 
and in this region shocks are not possible. They are, however, possible beyond r s = r* = r^^T" 1 , 
where the baryon Lorentz factor has saturated to the value T s = 77*. Any initial variations on a 
timescale t v will lead to shocks at a radius r s h given by equation (|i~3|) with T s ~ 77*, which will be 
located above r s = r#. There is thus a maximum possible internal shock radius for any given £„, 

r shM = 2.17 x l0 l ^rjlr o = 2.17 x 10 14 (L 52 Mir o y) 1/2 ^ cm , (15) 

which, unless £ v is large or the external density is very large, can still be smaller than the radius 
where external shocks are expected. The range where internal shocks can occur is shown as the 
horizontally or vertically striped region in Figure 1. 



-7- 



The internal shocks in the wind can dissipate a fraction of the terminal kinetic energy 
luminosity above the saturation radius r s . For a mechanical efficiency e s h = 10 _1 e s / l _i of 
conversion of kinetic energy into random energy which can be radiated, the shock luminosity 
in the radiative regime is 

L s h _ ( lQ~ l e sh -i , for r] < 77*, r > r s ; 

L 1 10 -1 e s / l _i(r7/^*) -1 , for rj > 77*, r > r s . 



4. Photosphere and Shock Spectra: Comptonization and Pair Formation 

The basic photospheric spectrum is that of a blackbody, xF x oc x s exp(— x/Q p h), with a 
thermal peak at 

x ph ~ 3Q ph < Xp ho = 30 o , (17) 

where Q p h,Q are given by equations (||]2]). Notice that at r s the level is Fpht = an d the 
spectral peak is at x p h ~ 30 o ~ 1, from equations @ @. For r p h > r s , both and Q p h 
decrease oc (r p h/r s )~ 2 / 3 . In an xF x or xL x diagram the thermal peak x p h (labeled with T in Figure 
2) moves down and to the left with a slope 1 as rj decreases. Before leaving the photosphere, 
however, the blackbody photons can act as seeds for scattering to higher energies, if there is a 
substantial amount of energy in scattering centers moving with characteristic speeds or energies 
larger than that of the emitting electrons (which are subrelativistic since Q' ph < Q' p = 0.03). 
Alfven waves generated by magnetic field reconnection or MHD turbulence can act as such centers 
(Thompson 1994). Alfven waves travel at speeds V w nearly the speed of light, with an equivalent 
comoving electron energy 6^ = kT' w lm e (? = (1/3) (7^, — 1) ~ (l/3)(I^/c 2 ) ^ (1/3), where 
V w /c ^ 1, and these waves can be efficiently damped for t s > 1. Alternatively, shocks which occur 
inside the photosphere may also induce Alfven waves. Repeated scattering on the Alfven waves 
acts in the same way as comptonization off hot electrons. The spectrum follows from conservation 
of photon number and conservation of energy. As seen in the observer frame, starting from seed 
photons at energy Q p h < O this yields a spectrum F x oc x° or xF x oc x, as the conserved photon 
number is scattered up in energy. The photons diffuse up and to the right on a slope oc x in the 
xF x diagram. From conservation of energy, the maximum energy they can reach is ~ 30 o ~ 1. 
However, if the energy in Alfven waves or turbulence is a fraction e w < 1 of the total L Q , the 
comptonized photosphere luminosity is 

Lphc = £u>L , (18) 

and the comptonized photosphere spectrum xF x oc x can extend only up to a break energy 

Xphc = min[e^30 o , {L phc / L ph )3Q ph ] < 30 o . (19) 

This is still much less than the wave equivalent energy in the lab- frame, 0^, = (kT^ / m e c 2 )ri ^ 
(l/3)r. Thus above this break energy, an increasingly smaller fraction of the total photons can be 
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scattered with spectrum xF x oc x up to a maximum energy in the lab frame 

x w ~ m w <, r . (20) 

As discussed by Thompson (1994) and in classical references on comptonization, such a spectrum 
xF x oc x° is naturally expected from the diffusion of photons out of bounded scattering regions 
(such as reconnection hot-spots or turbulent cells in this case). We show in Figure 2 the 
comptonized photosphere component (labeled PHC), assuming a turbulent wave energy level 
e w = 10 L G , for various values of r], while in Figure 3 a lower value e w = 10 -2 is assumed. 

Internal shocks outside the photosphere, expected in the coasting regime if the outflow is 
unsteady, provide a significant nonthermal component of the spectrum. The primary energy loss 
mechanism in the shocks is synchrotron radiation, or inverse Compton (IC). It is common to 
assume that at the shocks the magnetic field energy density is some fraction es of the equipartition 
value with the outflow. The dimensionless field at the base of the flow is then 

xbo = {B/B Q ) = (2L e B /r 2 c) 1 / 2 B Q l = ll)' 2 ^^ , (21) 

where Bq = lum^e" j eh = AAA x 10 13 G is the critical field. Thus the dimensionless comoving field 
at the shock is 

x'ssh = XBoirsh/ro)-^- 1 = l(T 7 r 2 3 L52 2 4 /2 CVr X • (22) 

Note that if the fields are turbulently generated and equipartition is with respect to L s h, instead 
of L Q , then from equation (|l6|), €b < min[e s /j, £sh(j)*/v)]- This is a time average value of B' 
over the duration of the shocks, and it neglects any time- varying compression factors associated 
with individual pulses. The observer-frame dimensionless synchrotron peak frequency in units of 
electron rest mass is x sy = (3/2)x^7^r, where the minimum electron random Lorentz factor in 
internal shocks 7 m ~ 0.9 x 10 3 e e is typically a fraction e e of the equipartition value 0.5m p /m e , 
remembering that internal shock collide at relative speeds T re i ~ 1. Thus the observed synchrotron 
peak is at 

x sy = x^T/r,*)- 2 = 1.26 x W^L^e^e 2 ^ 1 ^ 1 , (23) 

where x* = 1.26 x I0~ 1 e i g 2 e 2 fj, 1 1 ^' 2 (YT )~ l / 2 £ ! ~ 1 . The comoving synchrotron cooling time is 
t' sy = 2.5 x lO^V^- 1 = 2.5 x 10" 8 r|L 52 1 e B 1 e ( : 1 /ii^ s, using equation (H). For e B < 1 the 
IC mechanism can become important for the MeV radiation. However, for €b not too far below 
equipartition values, the IC losses occur on timescales comparable or longer than synchrotron, 
and produce photons well above the MeV range where breaks and anomalous slopes occur. The 
comoving expansion time is t' ex = r/cT = 0.65r2/ii£u s, and the ratio t' sy /t' ex exceeds unity only 
for rather large T <; 3 x 10 3 L^ 5 e^ 5 ej^ fi 1 1 ^ 5 ^ v l ^ if 7 > 7m- Thus the electrons are in the 
radiative regime above the synchrotron break x sy in all cases considered, and also for a large range 
below it. Above the break x sy the synchrotron spectrum is then xF x oc x°, or more generally 
xF x oc x~( p ~ 2 ^ 2 , from the power law electrons N^) oc 7 _p above 7 = 0.9 x 10 3 e e produced by 
Fermi acceleration in the shocks (for the rest of the discussion we assume p = 2 or oc x° above 
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the break as an example). Below x sy one has xF x oc x 1 / 2 , down to a synchrotron self-absorption 
frequency 

x a < 10^ 6 lT^b^ 6 1^^ /B , (24) 
where the equality applies if xF x oc x 1 / 2 down to x a , and the inequality applies if there is a 
transition to an adiabatic regime oc x 4//3 before reaching x a (for very high T or £ v ). Schematic 
synchrotron spectra are shown in Figure 2 from shocks where e s h = 10 _1 , and in Figure 3 (left 
panel) where e s h = 3 x 1CP 3 , assuming p = 2. These curves, labeled S, show the break at x sy and 
a radiative slope 1/2 below that down to x a . For cases where an adiabatic regime is achieved at 
frequencies below x sy but above x a , the synchrotron slope would steepen to 4/3 and x a moves 
further down in frequency. 

Pair breakdown via 77 — ► can occur when the comoving synchrotron luminosity from the 
internal shocks is large enough, which is the case over a non-negligible range of parameter space. 
The comoving compactness parameter is £' = 7iL(Tt(t /r), where = (aL s h/£irr 2 h m e c 3 r 2 ) is the 
comoving photon density and a is the fraction above threshold. For typical synchrotron spectra 
peaking at x sy and with xF x oc x° above that, a fraction above threshold in the comoving frame 
a = 0.3a_.3 is a typical value for the high £' cases. Thus 

£' = (aL sh a T /87rm e c 3 t v r 5 ) ~ 3 x ltfrj 5 L S2 £ a h-ia-.3^ 1 ■ (25) 

For values of T > r# ~ 3.1x ^{L^sh-i®-.?,^ 1 ^ 1 ) 1 ^ the compactness £' < 1 and pair formation 
does not occur. This corresponds to shock radii {r s h/r )i> > 0.7 x ^L^sh-i^-sf^i ^2 3 - Below 
that, in the range r s < r s h < r 8 hi>, one has £' ^ 1, and pair breakdown rapidly adds to the opacity 
and cooling. This leads to a self-regulating pair plasma whose scattering optical depth t s ~ 
several, with a characteristic mean e ± energy & c ~ m e c 2 /r s ~ 10 -1 (e.g. Svensson 1987, Ghisellini 
Sz Celotti 1999). Clearly, t s cannot be too large, otherwise advection and adiabatic cooling would 
dominate over comptonization and diffusion. At such nonrelativistic energies, cyclotron radiation 
produces seed photons at harmonics whose energy is much lower than x sy , but repeated scattering 
on the much hotter produces a comptonized spectrum F x > oc x'° up to the maximum energy 
Q' c . For r = r\ (and £' > 1), one then has an observer- frame characteristic break energy 

x c ~ io _1 r ~ ior 2 , (26) 

and below that a spectrum xF x oc x. The pair Comptonized luminosity is comparable to the 
(initial) synchrotron shock luminosity in the absence of pairs, L C /L Q = L s h/L given by equation 
(|l~6|) , since the self-regulation of the comptonizing pair plasma is achieved by a time-averaged 
balance between energy dissipation by the shocks and radiative losses. 

Above the break x c , one would expect a drop-off of the spectrum, in the absence of other 
effects. However, since internal shocks occur in the coasting regime and self-regulating pair 
breakdown tends to maintain a moderate scattering optical depth, reconnection and MHD 
turbulence may arise here too, leading to Alfven waves of much higher characteristic energies than 
that of the pairs, which could lead to a flatter power law spectrum xF x oc x° extending above x c 
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up to energies x w similar to that given in equation (^). In Figure 2 the top row shows cases where 
£' < 1 and hence only the photosphere (with thermal T and nonthermal PHC) and the shock 
synchrotron (S) components are present. The middle row shows the marginal cases where £' = 1 
and hence besides the photosphere T and PHC components one expects the synchrotron S and 
the pair breakdown comptonized C components to be comparably important. The bottom row of 
Figure 2 shows cases with £' ^> 1, where pair breakdown is so important as to completely replace 
the synchrotron component S with a self-regulated comptonizing pair plasma component C. Both 
the S and C components have a luminosity level given by a shock efficiency (e.g. Kumar 1999) 
e s h = 10 _1 in Figure 2, with the effects of a lower value e s h = 3 x 1CP 3 shown in Figure 3 (left). 

5. Discussion 

Within the framework of the standard internal shock model, we have analyzed the observable 
effects of the two major radiating regions, the photosphere and the internal shocks, which are 
expected to contribute to the flux from an unsteady fireball outflow model of GRB. Note that, 
as in most GRB models, the spherically symmetric assumption applies to any conical outflow 
whose opening half-angle is <; r . (We note also that, if the cone angle were very narrow, 
transverse pressure gradients would cause significant departures from radial outflow; under these 
circumstances, the r-dependences would change, although the qualitative features would not be 
substantially altered.) We have purposely left out of our discussion any radiation from an external 
shock, which is expected to occur at radii beyond those considered here (and which can add other 
radiation components, especially a long term afterglow). 

The standard internal shock model of GRB is generally assumed to produce its observed 
nonthermal radiation by the synchrotron (or possibly inverse Compton) process. Here, in addition 
to synchrotron we have also considered in more detail the role of the outflow photosphere and 
of possible nonthermal spectral distortions in it, as well as the role of pair breakdown in shocks 
with very high comoving luminosity. In a diagram (Figure 1) of radius vs. dimensionless entropy 
r\ = L/Mc 2 , the regions where the internal shocks are dominated by synchrotron radiation (and 
pair breakdown is unimportant, £' < 1) are shown by the vertically striped area S. The line where 
the compactness parameter £' = 1 runs parallel to the line r^ h for the photosphere in the coasting 
regime, in the same figure. The region where internal shocks are dominated by pair-breakdown, 
£' >, is given by the horizontally striped region C in Figure 1, where the shock spectrum is 
dominated by comptonizing pairs. To the left and below the photospheric lines, shocks would 
occur at high optical depths and their spectrum would be thermalized, adding to the purely 
thermal (T) and non-thermal (PHC) components emerging form the photosphere. These various 
spectral components are shown in Figures 2 and 3 in a power per decade xL x vs. x plot, where x 
is photon energy in electron rest mass units. 

In our earlier papers on fireball shock models of GRB and most subsequent related work, the 
role of photospheres and pair breakdown was briefly considered, but until recently the observations 
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did not appear to provide much support for their being important. The need for a non-thermal 
spectrum continues to be a strong argument for shocks and a synchrotron component, while for 
the less problematic case of large baryon loads the photospheres occur in the coasting regime, 
where their observer frame thermal luminosity is drastically weakened by adiabatic cooling. This 
is confirmed by the spectra of Figures 2 and 3, where for moderate to large baryon loads (low 
rj) and moderate variability > 10 2 or t v <; 0.3 s the thermal peaks T are strongly suppressed, 
especially in the more "conservative" region going part way above and below from the central 
and right-of-center panels of Figure 2. Pair breakdown is also a phenomenon considered in earlier 
papers (c.f. also Pilla & Loeb 1998, Papathanassiou & Meszaros 1996), which received less 
attention than it deserves because its importance appears to be restricted to a relatively narrow 
region of parameter space. This is illustrated in Figure 1, where one can compare the vertically 
striped region £' < 1 labeled S versus the narrowish horizontally striped region £' > 1 labeled C. 
This region would be even narrower if one normalized to L a ~ 10 50 erg s -1 as opposed to L52. 

However, fresh motivation for reconsidering the role of photospheres and pairs is provided 
by the evidence, in a non-negligible fraction of bursts, of low energy (1-10 keV) spectral slopes 
steeper than 1/3 in energy (or 4/3 in logxL-,;) and in some cases an X-ray excess above the power 
law extrapolation from higher energies, as well as the ubiquitousness of observed break energies 
clustering between 50-1000 keV (discussed in §1 and references there). A look at Figures 2 and 
3 shows that allowing a larger role to photospheres and comptonized pairs provides a way of 
addressing these observational trends. If photospheres are rather more important than shocks 
(e.g. due to some preference for very high 77, weakly varying outflows or low shock efficiencies) the 
thermal component can provide low energy slopes as steep as oc x 3 in xL x (while flatter slopes 
can be achieved through integration or distributions). The same can explain an X-ray excess 
well above the power law extrapolation from above. Reasonable break energies can be obtained 
from either the synchrotron or pair comptonization mechanisms in the shock, but they depend on 
r = 77, and unless the range of 77 is narrow they would not necessarily cluster between 50-1000 
keV. In the case of pair comptonization they also tend to be a bit high, unless the equilibrium pair 
temperature is @' c ^ 0.03. A preferred break could be attributed to a photospheric peak, provided 
baryon loads are very low, 77 > 77* ~ 10 3 in all cases; this still requires a strong shock synchrotron 
component, or possibly Alfven wave comptonization in the photosphere, to explain the high 
energy power law spectra. It would also imply a pronounced upward change of slope above the 
break, from the thermal peak to a flatter power law in all bursts where a break is observed, unless 
the shock or Alfven wave scattering always produces a luminosity comparable to the thermal 
photosphere. Preferred break energies arise naturally if Thompson's (1994) proposed mechanism 
of comptonization by Alfven wave damping in the photosphere is taken at face value and has good 
efficiency, leading to a source frame break around 0.5 MeV. Note that the comptonized spectral 
slopes discussed here (either from thermal pairs in shocks or Alfven waves in the photosphere) are 
nominally 1 and in xF x , but for simple time-dependent calculations an evolving slope is expected 
(e.g. Liang et al, 1999, Thompson 1994). In reality the actual time dependence for an unsteady 
outflow leading to shocks, pair breakdown and comptonization could be more complicated. 
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If both a photospheric and a shock component are detected, one would expect the thermal 
photospheric luminosity (and its non-thermal part, if present) to vary on similar timescales as 
the nonthermal synchrotron or pair comptonized shock component (unless the shock efficiency is 
radius dependent, or unless one or both are beyond r = r Q r] 2 , in which case r/crf imposes a lower 
limit on the corresponding variability timescale). But even if the bolometric luminosity varies on 
the same timescale, the luminosity in a given band (e.g. BATSE) probably varies differently, since 
the thermal peak energy is oc L 1 / 4 and falls off steeply on either side, while the synchrotron peak 
energy varies oc B'r\ oc L 3//2 and falls off more slowly. The pair comptonized break energy on the 
other hand varies as rj oc L, and also drops off slowly on the low energy side, or on both sides 
if scattering off waves is present in the high energy side. If observations at high time resolution 
become possible in X-ray or optical during the burst (as opposed to the afterglow), we would 
expect (c.f. Figure 2) the bursts with shorter time structure (low £„) to be more suppressed at 
these wavelengths, compared to those with longer variability timescales. 

In summary, there are several plausible mechanisms for producing a preferred energy break, 
which rely on internal properties of the outflow. In those bursts with low energy slopes steeper 
than implied by synchrotron, the prominence of the photosphere is a likely explanation, in 
which case its luminosity could in some cases vary differently from the higher energy power law 
component. Intrinsically high luminosity bursts, where pair breakdown is inferred, would be 
predicted to have generally harder high energy power slopes than in lower luminosity bursts where 
synchrotron provides the high energy slope. If photospheric comptonization on Alfven waves is 
responsible for the high energy power law slopes, the thermal peak and the power law should vary 
together in time. In this case a straightforward prediction is that the thermal peak photon energy 
and the intrinsic luminosity allow one to determine the expected maximum photon energy x w ~ T. 

This research has been supported by NASA NAG5-2857, NSF PHY94-07194 and the Royal 
Society. 
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Fig. 1. — Radius vs. rj = L/Mc 2 dia gram showing the location of the photosphere r p h, the 
saturation radius r s , and the shock radii r s h for various variability timescales £„. The regions 
where the shock spectrum is dominated by synchrotron (S) or comptonization in a pair plasma (C) 
correspond to a comoving compactness £' < 1 or £' > 1. 
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Fig. 2. — Luminosity per decade xL x vs. x = hu/m e c 2 for different values of 77 increasing upwards, 
and different variability timescales increasing to the right, using e w = e s h = 10 _1 . Pair formation 
is unimportant in the top row and important in the bottom row, being marginal in the middle 
row. T: thermal photosphere, PHC: photospheric comptonized component; S: shock synchrotron; 
C: shock pair dominated comptonized component. 
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Fig. 3. — Luminosity per decade xL x vs. x for two of the spectra in Figure 2 but with different 
shock synchrotron and pair comptonized components or lower comptonized photosphere component. 
Left panel: e w = lO" 2 ,^ = 3 x 1CT 3 . Right panel: e w = lO -2 ,^ = 1CT 1 . 



